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Abstract. Modeling most of the tests of general relativity requires to know the 
function relating light travel time to the coordinate time of reception and to the 
spatial coordinates of the emitter and the receiver. We call such a function the 
reception time transfer function. Of course, an emission time transfer function 
may as well be considered. We present here a recursive procedure enabling 
to expand each time transfer function into a perturbative series of ascending 
powers of the Newtonian gravitational constant G (general post-Minkowskian 
expansion). Our method is self-sufficient, in the sense that neither the integration 
of null geodesic equations nor the determination of Synge's world function are 
necessary. To illustrate the method, the time transfer function of a three- 
parameter family of static, spherically symmetric metrics is derived within the 
post-linear approximation. 
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1. Introduction 

In many tests of general relativity, it is crucial to explicitly know the function relating 
light travel time to the coordinate time of reception (or emission) and to the spatial 
coordinates of the emitter and the receiver. As the accuracy of measurements improves 
over the years, it will be soon indispensable to take into account the effects on the 
propagation of light predicted within the post-linear approximation of the relativistic 
theories of gravity. It will be the case, e.g., if experiments like the Laser Astrometric 
Test of Relativity (LATOR) [T] or the Astrodynamical Space Test of Relativity using 
Optical Devices (ASTROD) [2 are carried out. 

In most of the studies devoted to determining light travel time, the relativistic 
effects are calculated by integrating the null geodesic equations (see, e.g., Refs. [3] 
and [4]). General results have been obtained within the linear regime [3 [6l [7l [8]. 
However, the implementation of this method comes up against two problems. First, the 
integration of the null geodesic equations leads to heavy calculations in the post-post- 
Minkowskian approximation, even in the simple case of a static, spherically symmetric 
space-time [9j [TOj [11] [12] . Second, once the general solution is obtained, it is not easy 
to deduce from it the null geodesic path which joins an emitter and a receiver that 
have specified spatial positions. 

An alternative procedure enabling to overcome these two difficulties has been 
recently developed in [T3l[Ill[in]: the method consists in determining Synge's world 
function n{xA,XB) for any pair of points-events xa and xb [H], and then deducing 
light travel time from the equation Q{xa, xb) = which expresses that xa and xb are 
linked by a null geodesic path. This procedure works well, but presents an unpleasant 
drawback: once the general post-Minkowskian expansion of the world function is 
known, obtaining the corresponding expansion of light travel time still requires a lot 
of additional calculations (cf. Ref. |15j). It is the purpose of this paper to present 
a stand-alone method which totally avoids performing the calculation of the world 
function. 

The present work is organized as follows. In section 2 we summarize the notations 
and conventions used in the paper. In section 3 we recall the definition of the reception 
time transfer function and of the emission time transfer function, and we show that 
each of these functions obeys to a Hamilton- Jacobi-like partial differential equation. 
In section 4 we define the reception time delay function and the emission time delay 
function, and we show that each of these functions satisfies an integro-differential 
equation. Using this result, we derive in section 5 the general post-Minkowskian 
expansions of the time transfer functions. In section 6 we focus on the case of a 
stationary space-time. In section 7 we apply our method to the gravitational field 
of a static, spherically symmetric body treated within the second post-Minkowskian 
approximation. We present some concluding remarks in section 8. 

2. Notations 

In this work the signature of the Lorentzian metric g is (-1-,—,—,—). We suppose 
that space-time is covered by some global coordinate system (z^) = {x'^,x). We put 
x'^ = ct, c being the speed of light in a vacuum. Greek indices run from to 3 and 
latin indices run from 1 to 3. Einstein convention on repeated indices is used here for 
expressions like as well as for expressions like A^^B^ or a^Ci. Bold letters denote 
ordered triples. Given two triples a = {a^,a?,a^) — (a') and b = {b^,b'^,b^) = (6'), 
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we use a.b to denote a' 6*. The quantity \a\ stands for the ordinary Euclidean norm 
of a: \a\ — {Sija^a^^/^ . When it seems necessary for the sake of legibihty, a quantity 
f{x) is denoted by {f)^ or [f]x- The indices in parentheses characterize the order of 
perturbation. These indices are set up or down, depending on the convenience. 
G is the Newtonian gravitational constant. 



3. Time transfer functions 

We assume that space-time is globally regular with the topology of M x R'^, that is 
without horizon. We suppose in addition that the coordinate system is chosen in such 
a way that the curves of equation x = xa are timelike for any given xa, which means 
that d / dx^ is a timelike vector field, i.e. goo > everywhere. In agreement with these 
assumptions, we suppose that the past null cone at a given point xb = ictB,XB) 
intersects the world line x — xa at one and only one point xa — ictA,XA)- The 
difference is — tA is the (coordinate) travel time of a light ray connecting the emission 
point Xa and the reception point xb- This quantity may be considered either as a 
function of the instant of reception tB and of xa, xb, or as a function of the instant of 
emission tA and of xa and xb- So it is possible to define two time transfer functions, 
%. and Tf. by putting 

tB-tA^ Tr{xA,tB,XB) Te{tA,XA,XB). (1) 

We call % the reception time transfer function and % the emission time transfer 
function. These functions are distinct except in a stationary space-time in which the 
coordinate system is chosen so that the metric does not depend on x". 

Let us denote by k'^ the vector dx^^/dC, tangent to the null geodesic path 
connecting xa and XB^ C being an arbitrary affine parameter of this geodesic. It 
was shown in [M] that the covariant components of fc'^ satisfy the relations 

ki\ _ dTr{xA,tB,XB) , . 

ko).-' dx^A ^ ^ 



h\ _ _ dTe{tA,XA,XB) , s 



and 



The following theorem can be inferred from equations Q and 

Theorem 1. The time transfer functions %.(xA,tB,XB) and Te(tA,XA,XB) 
satisfy the Hamilton- J acobi-like equations 

d% 



- c%, xa) + 2c5*(4 - cT,, xa) „ , 

dx\ 

d% dTr 



c'f^(a:'^-c%,XA)-^-f^O (4) 
ox A dx\ 



dT 

{x\ + c%,XB)-2cg°\x''A+cT„XB) " 



dxB 



+ c2 5'-' (xUcre, 03^)1^^=0, (5) 
c>x'b dx^g 

X Note that the order of the arguments of Tr in equation JTJ is different from the one given in Refs. 
[Til and [T5l . 

§ Equations (|4]l and l|5} can also be derived from the general-relativistic version of Fermat's principle 
established in Ref. |17| . 
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respectively. 

Proof of theorem 1. The covariant components of the vector tangent to Tab 
at XA satisfy the equation 

{g'-^Wx, = 0. (6) 

Dividing equation (iniside by side by [(/co)^]^, and then taking equation ^ into 
account yield equation (HlljJ. The same reasoning using {g'^'^k^k^)^^ = and equation 
([3]) leads to equation ([5]). 

In what follows, we give detailed proofs only for the reception time transfer 
function because similar procedures may be applied to the emission time transfer 
function. 



4. Integro-differential equations satisfied by the time delay functions 

Henceforth we suppose that the metric takes the form 

gfj.,^ = Vfiiy + h^iv (7) 
throughout space-time, with rj^i, = diag(l, — 1, — 1, — 1). Then the contravariant 
components of the metric may be decomposed as 

g^-'^ ^ r1^-'^ k^"' , (8) 
the quantities fc^'' being determined by the relations 

fcM- + T^^Pry'^-Zip, + fj^'Php^k'"' ^ 0. (9) 
According to equation ([7]) the reception time transfer function may be written as 

%{XA, tB,XB) = -RaB + -Ar{xA, tB,XB), (10) 

c c 

where 

Rab = \xb - xa\ (11) 

and Aj,{xa, ts, xg) is of the order of the gravitational perturbation h^^i,. The function 
lS.j./c defined by equation (|10p may be called the reception time delay function. It is 
well known that this quantity is > in Schwarzschild space-time, which explains its 
designation. 

Replace now xa by a variable x and consider f^, xb as fixed parameters. Inserting 
Tr{x, tB, Xb) = \xb — x\/c + Ar{x, tB, X b) / c mto equation (j4]) taken at x instead of 
Xa, we get an equation which may be written in the form 

2N^'-^^:^^0^^-Wix,tB,XB), (12) 



where iV* is defined as 



X 



B 



N' = 

\xb - x\ 

and W{x,tB,XB) is given by 

W{x,tB,XB) = k^^{x^) - 2N'k°\x^) + N'N^k'\x-) + 2 [k°\x^) ~ Wk'\x^)\ 

dAr-{x,tB,XB) , r J, , si fdArdAr' 



X 



dx* '■ ■' \ dx^ dx^ 



(x,tB,XB) 



II The covariant component feo of a null vector k cannot vanish in the chosen coordinate system since 
the vector d/dx'^ is assumed to be timelike everywhere. 
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X- being the point-event defined by 

= {x% - \xb - x\- Ar{x, tB, xb),x). (14) 

Since a; is a free variable, consider the case where x is varying along the straight 
segment joining Xj^ and Xb- Then we have 

= N\b, (15) 
where is by definition 

N\b = (16) 

and we can put 

x = z-{X), (17) 

where 

z.{X)=xb-XRabNab, 0<A<1. (18) 

A straightforward calculation shows that the total derivative of Ar{z^ (A), tB , xb) 
with respect to (w.r.t.) A is given by 

— Ar{z^iX),tB,XB) = -RABN\B-g^iZ-W,tB:XB), (19) 

where dAr{z-{X),tB,XB)/dx'^ denotes the partial derivative of Ar{x,tB,XB) w.r.t. 
X* taken at a; = Z-{X). 

Taking equation (jlSp into account, and then comparing equation (|19p with 
equation (H^), it may be seen that Ar(2;_ (A), t^, £c_b) is governed by the differential 
equation 

^Ar{z_iX),tB,XB) = ^RABWiz_iX),tB,XB) (20) 

with the boundary condition 

Ariz^{0),tB,XB)=0 (21) 

which follows from the obvious requirement Ar{xB, tB,XB) — and from xb — Z-{0). 
As a consequence Ar{z^{X),tB,XB) is such that 



(22) 



Ar{z^{X),tB,XB)^^RAB I W {z_{X') , tg , X B)dX' . 

^ Jo 

Noting that ^-(1) = xa, we deduce a theorem as follows from Eqs. ([22|) . (fT3 
and dl 



Theorem 2. The reception time- delay function Ar defined by equation ilO\ 
satisfies the integro- differential equation 



Ar{xA, tB,XB) - ^RaB ^ - ^N^Bk"' + A^AB^^isfc'') „ 



{k--Ni,k^^) J^iz^iX),tB,XB) 



[r^^^ +k^^{z^X))] 



z_(A) dx 

dAr dAr 



dX, (23) 



(z_(A),i_B,X_B) 



dx'^ dx^ 

where z-iX) is the point-event defined by 

z_(A) = ix% - XRab - Ariz_iX),tB,XB),z.iX)), (24) 
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Z-{X) being given by equation hl8\) . 

Of course, if we define the emission time delay fmiction Ae(iA; ^c^, a;B)/c by the 
equation 

Te{tA,XA,XB) = -RaB + -/^e{tA,XA,XB). (25) 

c c 
a similar theorem may be stated for Ag. 

Theorem 3. The emission time-delay function Ae defined by equation I125\) 
satisfies the integro- differential equation 



1 r 

A,itA,XA, xb) = -Rab (fc"" - 2N\Bk°' + N\sN'^gk'=y^ 

-2{kO^-NUk^^)^^J^itA,XA,.,i,)) 



dfi, (26) 



where z+(/i) is the point- event defined by 

Z+{fl) = {xa + fJ.RAB + Ae{tA, XA, Z+{p)), (27) 

with 

z+{fi) = XA + uRabNab, (28) 

and dAe{tA,XA, z+{^))/dx'' denotes the partial derivative of Ae{tA,XA,x) w.r.t. x* 
taken at x ^ 2;_|_(/i). 

Integro-differential equations ((23|) and (|26p are very convenient to obtain the 
general post-Minkowskian expansion of each of the time transfer functions, as we shall 
see in the next section. 

5. General post-Minkowskian expansion of time delay functions 

Henceforth we suppose that each perturbation term /i^^ is represented at any point x 
by a series in ascending powers of the Newtonian gravitational constant G 

oo 

V(a;,G) = ^G"g(";(x). (29) 
The concomitant expansion of fc^" is then given by 

oo 

fc^"^(x,G) = 5]G"gf„yx), (30) 
where the set of quantities g'f'^, can be recursively determined by using the relations 



= -v^'v-'^g^^J (31) 



and 

n-l 

J(n) ~ 'I 'I 'I ypcr y(n-p) i-^^J 



p=l 

for n> 2. 
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As a consequence, the reception time delay function admits the expansion 

oo 

Arix, tB,XB,G) = J2 G"A(")(a;, tg^XB). (33) 

n=l 

It follows from equations ([50)1 and ([M)) that each term k^'^ involved in the right 
hand side of equation (j23[) may be written as 

oo 

F''(F_(A),G) = ^G"5[:r)(4 - ^Rab " A,(^_(A),ts,XB,G),^_(A)). (34) 

n=l 

The general post-Minkowskian expansion of A;'"'(zL (A), G) is obtained by 
substituting for from equation (I33p into equation (I34p . and then performing 
Taylor's expansion about the point Z-(A) defined by 

z.iX)^ix%-XRAB,z_iX)). (35) 

Let us put 

k 

l[Ai'^ + 'Hx,tB,XB) 

where Zi, /2, ^fc are either positive integers or zero (m > 1 and 1 < fc < to). A 
straightforward calculation yields 



<flj^^''Ux,tB,Xi 



k\ 



E 

iH hlk=m — k 



(36) 



fc^"^ (z_ ( A) , G) = £ G"g':'^„) (^_ ( A) , , a^B) , 

n=l 

where the quantities 5^'^^j(z_(A), t^, cc^) are given by 
?^!'^l)(^-(A),^B,a;B) =g[;'^(2;-(A)) 

and 

9-\n) (A) , is , a^s) = .glt^) (^- (A)) 



9' 



k„k^^ 

■0\k 



{dx' 



(37) 

(38) 
(39) 



J2_(A) 



+ ^^$('».fc)(z_(A),iB,a;B) 

m=:l fc— 1 

for n > 2. Substituting for fc^'^(z_(A), G) from equation ([57]) into equation we 
get the theorem which follows. 

Theorem 4. When the perturbation term h^^i, in the metric is represented by 
expansion i29\) . the reception time-delay function A^ is given by equation I133\) . namely 

oo 

Ar{x, tB,XB,G) = G"A(")(a;, tB,XB), 



where 

Al.^\xA,tB,XB) = ^RaB 
Al^\xA,tB,XB) = ^RaB 

+ 2 



{z^{X),tB,XB) 

dX 



dX, 



(40) 



?(i) ^^ab9{\) 
aA« 9A«' 



9a;* 9a;^ 



(z_(A),tB.a;B) 



(41) 



Time transfer functions 



B,Xb) 



n-l 



(2-(A),tf3,a;ij) 

(2;_(A),ts,a;s) 



E 



77- 



n-2 

+ '^g-(p-,iz~W,tB,XB) 



- (z-{X),tB ,xb) 



9=1 



aA}."^ 9A, 



dx"^ dx^ 



dX 



(42) 



for n > 3, t/ie quantities dl'^n'j being defined by equations i38\) and i39\) . 

It may be noted that the integral expressions occurring in Eqs. (|40I) - P^ are hne 
integrals taken along the zeroth-order null geodesic of parametric equation x — ^-(A), 
where z_(A) is defined by equation (|55|) . 



A similar reasoning works for the emission time transfer function, which admits 
the expansion 



AeitA, XA, X,G)^Y1 G''A^J'HtA,XA, x). 



(43) 



We define the functions ^e"^'''\tA, xa, x) as 



llA(^'+'\tA,XA,x) 



^i"^'''HtA,XA,x)^^^ ^ 

h-i i-lk=m-k 

where li,l2,---,lk are either positive integers or zero. Then, defining ^+(/i) as the 
point-event 

Z+{^l) = {x°^ + nRAB,Z+{fi)), 

where z+(/i) is given by equation (I28p . we put 

9t1i)i.tA,XA,z+{fi)) = ,g['j'^(z+(^)) 

and 

d+\n)itA,XA,Z+{^i)) ^ .g[',^^(z+(^)) 

r ^, 



, (44) 



(45) 
(46) 



n — 1 m 



+ EE'^'e'"''^(*-4,a:A,^+(M)) 

m— 1 k—1 

for n > 2. Thus we can formulate a theorem as foUows. 



gk 



(dx' 



0\k 



(47) 



J z+(At) 
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Theorem 5. On the asumption of Theorem 4, the emission time transfer function 



IS 



given by equation ^3\ l, where 



A';^'>{tA,XA,XB) 



Rab 
Rab 
-2 



1 ^ 



(48) 



(tA,a;A,z+(Ai)) 



ri- 



all) - Kb9^ 



-{tA,XA,Z+{^l)) 



- {tA,XA,Z+(tJI.)) 



(49) 



and 

A^J'\tA,XA,XB) 



— t^Rab 



9+{n 



ri-1 



{tA,XA,Z+{^l)) dx^ 



-{tA,XA,Z+{fJ.)) 



n-1 

p=i 

n-2 



(P) ;3a("^p) 



dA)F> dA 



dx'- dx^ 



n—p—l 

J29+{p)i*A,XA,Z+{ll)) 
p=l q=l 



- (tA,XA.Z-i-{fi)} 



dx^ 



dx^ 



dn (50) 



(tA,a;A,z+(A')) 

for n > 3, the quantities g'^(^n) ^^''^5 defined by equations and ^42^- 

The integrals are now taken along the zeroth-order null geodesic defined by the 
parametric equation x — z^{fi). 

As far as we know, the above results are new for n > 2. On the other hand, 
the expressions obtained here for Af^' and Ai""^' are equivalent to the well-known 
expressions derived by other methods. Taking equation pip into account, it may be 
seen that the formulae (|40| and (|48|) coincide with the expressions obtained, e.g., in 
Ref. [14] for cT^^^"* and cTe^\ respectively. 



6. Case of a stationary space-time 

In the case of a stationary space-time, one can choose coordinates (a;") such that the 
components of the metric do not depend on a:°. Then the two time transfer functions 
reduce to a single one, so that we can write: 

%{xA,tB,XB) = X,{tA,XA,XB) = T{xa,Xb)- (51) 

As a consequence, the functions A^ and Ag are identical and depend only on xa 
and xb, which implies that 

Ai"^{xA,tB.XB) = Ai"'>{tA,XA,XB) ^ A^''\xa,xb) (52) 

for any n > 1. 
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Moreover, the stationary character of the metric imphes that 

r_\^^{z_{\),tB,XB)^9'^:^{z.{x)) 

and 

So equations (HDl) and P5)) - ([50)) simphfy. For example, equations 

(|49p may be written as 

»i 



and 

A^'Hxa,xb) 
-2 



1 



-Rab 



:Rab 



gfl) - 2N\Bgtl) + N'ABN'AB9n) , , 
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(53) 

(54) 
and 

(55) 



5(1) ~ ^^ABff(l) 



(a3A:Z+(M)) 



(56) 



respectively. 

A straightforward calculation shows that making the change of variable 
A = — /i 

transforms the right hand side of equations ([55]) and ([55)1 into the expressions deduced 
from equations (|^(])) and ([H]) , respectively when equations ([5^ and ([55]) are taken into 
account. It is worthy of note that the integrals reduce in this case to line integrals 
taken along the segment joining xa and xb- 



7. Application to a static, spherically symmetric space-time 

To illustrate the previous results, let us consider the gravitational field outside a static, 
spherically symmetric body of mass M . Choosing spatial quasi Cartesian isotropic 
coordinates and putting r — |a;|, we suppose that the metric has the form 

^ 2 A 2GM ,^G2m2 \ , 



GM S^G^AP 
2' 



l + 2-f— + -^-J^ + ...]S^Jdx'dx', (57) 



where (3 and 7 are the usual post-Newtonian parameters, ^ is a post-post-Newtonian 
parameter and -t-... means that terms of order G^ are neglected (/3 = 7 = 5 = 1 in 
general relativity). Furthermore, we suppose that points xa and xb are such that the 
geodesic path connecting them is entirely outside the body. We use the notations 

rA = \xA\, rB = \xB\. 
The contravariant components of the metric are given by 

9^.^ = '^^ 9"ih=0^ 9l-,=^lf-/', (58) 

5^=2(2-/?)^, 5?^)=0, 9l) = q5-^^^)^,5-. (59) 
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Taking equations (|58p into account, we immediately deduce from equation (j55p 
that A(i) is given by the weU-known formula 

which is equivalent to the expression of the time delay found by Shapiro [18j . 
Substituting now equations ([58| and (|59|) into equation ((56|) . we get 



(61) 

We deduce from equation ([501) that 

A.i,^(..,.,(,))^(. + 1)^^ (62) 



and 



rA|2;+(/i)| + hxa-{xb - xa) + r\ 

(63) 



Substituting equations and (|63p into equation (pT|) . and then noting that 
d/i arccos(n^.nB) 



rArsy/l^^inZnB^^' 



(64) 



where 



riA = — , ns = , (65) 

rA TB 



we obtain by a straightforward calculation 

(8 - 4/? + 87 + 3(5) arccos(nA.nB) (1 + 7)^ 



A^^^(a;A,a;B) = 



rATB 



4^1 - {nA-riBY 1 + (ma-"b) 



(66) 

Finally, using the formulae ((60)) and ((66)) . we get an expression as follows for the 
time transfer function 

^, . Rab , (7 + 1)0^/ (rA + TB+RAB 

T{xa,xb) = + 5 In — 

c \rA + rB-RAB, 

S-4/3 + 87 + 3(5)arccos(nA.nB) (1+7)^ 



^ G-'M^ Rab 
& rArB 



4^/1 - (nyi-Wfi)^ l + (nA.nB) 

+ 0(G^). (67) 

We recover a result previously derived by different approaches [TOl E] (see also 
Refs. [IT] and [H] in the case where /3 = 7 = (5 = 1). 
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8. Concluding remarks 

Equations ([40 |) - (|42|) and ([48|l - (|50l) are the main results of this paper. These 
equations show that the time transfer functions can be obtained within the nth 
post-Minkowskian approximation by a recursive procedure which spares the trouble 
of solving the geodesic differential equations and avoids determining Synge's world 
function. It is remarkable that any nth-order perturbation term is an integral taken 
along a zeroth-order null straight line. The derivation of the time transfer function 
performed here for a static, spherically symmetric space-time is significantly more 
simple than the calculation carried out in Ref. |14j . 

As a final remark, it may be noted that since the time transfer functions are 
sufficient to determine the deflection and the frequency shift of a light signal, the 
systematic method developed here will be very convenient to tackle the relativistic 
problems raised by highly accurate astrometry and time/frequency metrology. 
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